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Abstract 

We have used an atomistic ah initio approach with no adjustable parameters to compute the 
lattice thermal conductivity of Sio.sGeo.s with a low concentration of embedded Si or Ge nanopar- 
ticles of diameters up to 4.4 nm. Through exact Green's function calculation of the nanoparticle 
scattering rates, we find that embedding Ge nanoparticles in Sio.sGeo.s provides 20% lower ther- 
mal conductivities than embedding Si nanoparticles. This contrasts with the Born approximation 
which predicts an equal amount of reduction for the two cases, irrespective of the sign of the mass 
difference. Despite these differences, we find that the Born approximation still performs remark- 
ably well, and it permits investigation of larger nanoparticle sizes, up to 60 nm in diameter, not 
feasible with the exact approach. 

PACS numbers: 66.70.Lm, 63.20.dk, 63.50.Gh 
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I. INTRODUCTION 



Nanoparticle Embedded in Alloy Thermoelectric (NEAT) materials have been proposed 
as a means of improving the thermoelectric properties of solid solutions. In a clear experi- 
mental demonstration of this concept, Kim et al. found a remarkable reduction of the ther- 
mal conductivity (k) of InGaAs upon the introduction of lattice matched ErAs nanoparticles, 
below the alloy limit, without any decrease of the thermoelectric power factor-. In some 
cases nanoparticles may also play an active role in increasing the power factor-. Naturally 
forming nanoinclusions are also thought to be at the core of the rather high thermoelectric 
figure of merit (ZT) of LAST (Lead-Antimony-Silver- Tellurium) materials^. Theoretically, 
the introduction of nanophases inside alloys to reduce k and improve ZT has been investi- 
gated using various approaches^"-. However, these approaches approximate wave scattering 
via a continuum description, and they rely on adjustable parameters. Our aim is to over- 
come these drawbacks by performing a parameter free atomistic calculation, which includes 
the nanoparticle scattering to all orders. We will show that the latter has a large influence 
on the results, which become highly asymmetric with respect to the sign of the scatterer's 
mass difference. (First order perturbation theory yields a quadratic, symmetric dependence.) 
Thus the question of whether it is better to embed lighter or heavier nanoparticles becomes 
very relevant in the light of the full order calculation. We provide the answer in the specific 
case of SiGe alloys, and we discuss its implications when developing novel NEAT materials. 

Our parameter free approach to compute thermal conductivity is based on an exact nu- 
merical solution of the linearized Boltzmann Transport Equation (BTE) for phonons^i"— . 
As a result of the computational challenge, many works still circumvent this solution by 
resorting to a host of approximations, most notably the relaxation time approximation^^. A 
few years ago some of us showed that it is possible to predict the lattice thermal conductivity 
of group IV single crystal semiconductors from first principles, i.e. using the fundamental 
physical constants as the sole experimental inputs^*^. In the present paper we extend our 
approach to the case of disordered solid solutions like Si2,.Gei_a;, including also embedded 
nanoparticles. This requires the ab initio calculation of elastic scattering rates due to both 
alloy disorder, and due to the nanoparticles. Here we use atomistic Green's function tech- 
niques to compute those rates to all orders, beyond the Born approximation. We will show 
that the full result may deviate noticeably from the Born approximation. 
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II. THEORY 



The thermal conductivity of a bulk material can be calculated as 

1 K, 



J2Mno + l)K\'hW,T,, (1) 

A 

where Vuc is the unit cell volume, the summation sign is a shorthand for = J^^ dq' 
(the integral is performed over the volume of the Brillouin zone), and A stands for the 
phonon branch index and wavevector, uj\ and are, respectively, the frequency 

and the group velocity along the z-direction of the corresponding phonon and is their 
occupation number. The t\ are scattering times that contain all the information about the 
non equilibrium phonon distribution. Details on the BTE and its solution have been given 
in^"— . We will just summarize it briefly. The equation to solve is 

r, = rl + r°A,, (2) 

where 

+ 

^A = X^r+^,^„(^AA"T-A" - Ga'T-A') + 
A'A" 

2^^-^'^"''^^^"'^^" ^ ^^^''^^') + Taa'^a'^a', (3) 

A'A" A' 

where, the Aa term takes into account coupling of non-equilibrium A phonon modes to other 
phonon modes (A', A") based on energy and momentum conservation, and 

l/rA° = E r+'A" + E ^^A-A'A" + E ^AA'. (4) 
A'A" A'A" A' 

The (+) and (-) symbols over the sums in Eq. H] indicate sums over A', A" for the two types 
of three phonon processes available, A ± A' ^ A". The meanings of the three-phonon terms. 



-'^AA'A"' ^^'^ ^AA" are given in Ref. lll|. The t\ were obtained through iterative solution of the 
BTE. For simplicity, the last term in Aa (Eq. |3]) has been neglected. This term vanishes 
for the nearly isotropic elastic scattering of low frequency phonons, which dominate the 
thermal conductivity. The effect of elastic scattering from nanoparticles and alloy disorder 
is included in the sum of Faa' in Eq. HI 

The main effect of a Si or Ge impurity or nanoparticle on lattice vibrations is through 
its mass difference compared to that of the host lattice. The dynamical equation for the 



displacements Ui is u;^Mu = /Cu, where K,ij = q^q^, is the interatomic force constants 
matrix, and My = Mi6ij is the mass diagonal matrix. When substitutional impurities of 
different mass are inserted, a diagonal matrix perturbation is added to the equation as 

V = — (M' — M)a;^ where M' represents mass matrix of the impure system. This pertur- 
bation is non zero only on the degrees of freedom associated with the scatterer. In the 
practical solution of the problem we work with the mass normalized matrices, K = M"^/C, 

V = — M^^V. It has been shown that the other perturbation term, corresponding to the 
differences in force constants, /C' — /C, has a much lesser effect on the thermal conductivity 
of SiGe alloys (about 10% of the totaU^.) In the case of nanoparticles, an estimation of 
its order of magnitude can be easily made. Rayleigh scattering due to mass difference or 
to differences in the elastic constants have similar expressions, except for the prefactor, 5*^. 
In the first case, this prefactor goes as ^mass ~ (^n. — Mm) /Mm- In the second case it 
is SqI ~ I " ^i"' 1 , where sub indexes m and n stand for matrix and nanoparticle re- 
spectively, and V is the speed of sound. The presence of the term is needed to ensure 
that SqI is zero if only the atomic masses change, but not the IPC's. For Sip.sGeo.s we have 
'5'mass ~ 0-4' whereas ~ 0.01 for L acoustic modes, and S^^^ ~ 0.0004 for transverse 
ones. Thus, the effect of different elastic constants is much smaller than that due to mass 
differences. Strain effects on the nanoparticle due to lattice mismatch could also induce ad- 
ditional scattering. In such case, a factor of order ~ 7(an — am)/am needs to be added to S^i, 
where 7 ~ 1 is the Griineisen constant, and a are the lattice constants of the nanoparticle 
and matrix materials respectively^. Addition of this term does not change the order of mag- 
nitude of Sgj, which remains smaller than ~ 0.02 and can be considered negligible compared 
with the mass difference effect. Taking into account the strain and IFC difference effects 
in an ab initio calculation would require the self consistent atomic relaxation of extremely 
large supercells, several times the size of the nanoparticles considered here. Therefore, given 
these difficulties, and the minor resulting effect on the total thermal conductivity, the IFC 
differences will be neglected here. 

The exact elastic scattering amplitudes due to a random distribution of independent 
scatterers in a homogeneous medium is Faa' = J2p f^^lx' with 

where is the volume fraction of scatterers of type p, Vp is the scatterer's volume, fl is the 
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volume into which the phonon eigenstates |A) are normahzed, and Tp(a;^) is the T matrix 
associated with the scatterer of type p^. We have adopted a virtual crystal approximation 
(VGA) model for the medium, where the interatomic force constants and atomic masses of 
pure Si and Ge crystals are averaged according to their relative concentrations in the alloy. 
The total alloy scattering for bulk Si2;Gei„2; is given by the concentration weighted sum of 
the scattering probabilities of a Si impurity in the VGA medium, and a Ge impurity in the 
VGA medium: Vf^^ = xTfy + (1 - a;)f^^,. 

The matrix T is defined in terms of the perturbation matrix V and the perturbed Green's 
function G"*" as T = V+ VG^V, which, after some algebraic manipulations using the orthog- 
onality and completeness of the eigenstates, can be expressed in terms of the unperturbed 
Green's function g"'"(a;^) as 

T{co') = [I-Yg^ico'r'V. (6) 
The integral form of the unperturbed Green's function g+(u;^) is given by 

mm^ (7) 

A 

where |A) are the eigenstates of the infinite unperturbed lattice and \i) is a local displacement 
of the ith degree of freedom in the direct lattice. For the numerical computation of g"'"(a;^) 
in Eq. [Tj we have employed the tetrahedron approach of Lambin and Vigneron^^. The total 
scattering rate due to the nanoparticles, appearing as the third term on the right hand side 
of Eq. m is efficiently computed using the optical theorem^^: 

l/r7 = ETaa' = ^2^H{m-^'m (8) 

A 

Most often in the literature, where the Born approximation is employed, the T matrix is 
replaced by the perturbation matrix V. This is justified by the expansion T ^ V + VgV + ..., 
valid for small perturbations. In this approximation one obtains 

1 ^ 2 n5 



itiiu"^) = lim 



= $]|5^e^ef |VA<,p5(a;A'-c.A) (9) 

A' k' 

Here, g2 = /p(l — M'/M^, Df is like a phonon density of states but weighted by the 
structure factor for the nanoparticle: = XliGnp being the number of unit 



cells making up the nanoparticle. The /' sum is only over those unit cells of the virtual 
crystal occupied by the atoms of the nanoparticle and Aq = q — q'. The are phonon 
eigenvectors for the k th atom in a unit cell. In the limit of a single atom impurity, Eq. M 
correctly reduces to the form derived previously by Tamura for the scattering rate of isotope 
impurities in cubic crystals^^. 

A consequence of the Born approximation is that the sign of the perturbation does not 
matter: a given percent of either increase or decrease of the scatterer's mass density with 
respect to the host's should produce the same result. This is not true when the exact T 
matrix is employed, and large differences can occur with respect to the Born approximation 
for large mass difference, as our results show below. 

III. RESULTS AND DISCUSSION 

We have first computed the thermal conductivity of Sio.sGeo.s at 300K and 800K. The 
perturbative approach for isotopic impurities^, where g2 ~ 10""^, has for decades also been 
used for alloys^, where one might question its validity since g2 in the alloy is several orders 
of magnitude larger (for Sio.sGco.s, g2 = 0.2). We have compared the Sio.sGeo.s scattering 
rates from Eq. (Q with those obtained using the full T-matrix method, and we find these 
to be close, as are the alloy thermal conductivities: h^qok = 10.62 (T-matrix) vs. 10.27 
(Born) W/m-K and ksook = 6.07 (T-matrix) vs. 6.0 (Born) W/m-K. Note that the alloy 
thermal conductivities are far lower than those of either bulk Si or bulk Ge because the alloy 
scattering is much stronger than the three- phonon scattering. As a result, A a in Eqs. ([2]) 
and ([3]) is small and tx ^ r°. These values are about 30% larger than the experimental 
values at the same concentration and temperatures (~ 7.5 — 8W/m-K at room temperature, 
~ 4.5 — 5 W/m-K at SOOK)^'^^'^^. There are several reasons for this. First, as already 
discussed at length in section II, our neglect of differences in force constants may lead to 
somewhat higher values. Also, experimental samples contain a certain amount of impurities 
and defects which also lower the thermal conductivity. In addition, there is a considerable 
spread in experimental results from different sources, which further attests to the various 
unknown factors present in experimental measurements of alloy samples. Finally, a recent 
first principles calculation of SiGe alloy thermal conductivity showed that the virtual crystal 
approach slightly underestimates phonon scattering in alloys^. 
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FIG. 1: Comparison of the scattering rates I/t^^^ normalized by nanoparticle volume fraction 
fp, due to Ge and Si nano particles of different size, in a Sio.sGeo.s matrix. The curves shown 
correspond to the LA phonon branch along direction (100). 

Our calculation shows that for single atom scatterers, low frequency phonons are well de- 
scribed by the Born approximation, with almost no difference between Si and Ge impurities. 
This is clearly seen in figure 1, for the case of the longitudinal acoustic branch. Although, 
the scattering rate of Ge and Si single atoms differ importantly at high frequency, however 
this does not lead to much difference in k, since k is dominated by low frequency phonons. 
It is only when we consider larger nanoparticles that differences become appreciable in the 
thermal conductivity. For nanoparticles of diameter 1.1 nm, containing 38 atoms, scattering 
rates differ considerably between the Si and Ge cases, already for u above 5 THz. This 
leads to a difference between the thermal conductivities of the corresponding composites. 
The difference becomes even more appreciable for larger nanoparticles (see Fig. 1). Fig. [2] 
compares the Born approximation result, with the exact scattering rates 1/t^^ for both a 
Si nano particle in Sio.sGco.s and a Ge nano particle in Sio.sGeo.s for a diameter of 3.3 nm, 
when the incident phonon direction is (100). 

At high frequency both the Si and Ge cases deviate considerably from the Born approxi- 
mation result (also shown in Fig. [2]). This occurs when the wavelengths become comparable 
to the size of the scatterers, so we are no longer in the Rayleigh regime. An earlier interpo- 
lation formula had been proposed to link between the Born and geometric regime scattering 
cross sections^iSSr^: 1/cr ^ 1/cgoom + 1/^Born- Fig. [3] and |4] show n versus nanoparticle size 
obtained using this approximated interpolation formula. There is an optimal nanoparticle 
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FIG. 2: Scattering rate due to Ge and Si nano particle of diameter 3.3 nm, in a Sio.5Geo.5 

matrix, normalized by nanoparticle volume fraction fp. The black dotted line is the Born ap- 
proximation result from Eq.(l9]). The curves correspond to the LA phonon branch along direction 
(100). 



size that minimizes thermal conductivity at a given concentration. These ab-initio curves 
confirm the simpler model predictions in Ref. Q, yielding a minimum for a diameter of a few 
nm, and a slow increase after that. 

The interpolated expression still makes use of the Born approximation, so it does not 
inform us of possible differences between heavier and lighter scatterers. We have compared 
those results with the ones obtained using the T-matrix computed scattering rates. The plot 
shows quantitative differences, but the trends are the same. The T-matrix approach is very 
computationally demanding: a 4.4 nm diameter particle, containing 2122 atoms is already 
at the limit of our computing capability. Therefore, we cannot assess the exact position of 
the minimum for the Ge or Si nanoparticle cases, although the graphs suggest that it may 
take place at a diameter between 5-10 nm. 

At equal nanoparticle size and concentration, the calculated n is always smaller for 
Ge (heavier) than for Si (lighter) nanoparticles. Their density difference with respect to 
Sio.sGeo.s is the same except for the sign, so in the framework of the Born approximation 
they should show an identical effect. This is clearly not the shown in Fig. |5l This 

figure shows the ratio between the conductivities of the two function of nanopar- 

ticle diameter: the thermal conductivity of the Ge nanoparticle case can be up to 20% lower 
than that of the Si nanoparticle case. The difference between the two cases highlights the 
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FIG. 3: Thermal conductivity vs. nanoparticle diameter at 1% (top), 3.4% (middle), and 5% 
(bottom) nanoparticle concentrations, for temperature 300K. Triangles: Born+geometrical inter- 
polation. Squares: T matrix calculation for Si nanoparticle. Circles: T matrix calculation for Ge 
nanoparticle. 

very different densities of the particles and the medium, close to 20%. The fact that Ge 
nanoparticles affect k more than Si nanoparticles is directly linked to their higher scattering 
rate at low frequency, visible in Fig. |2j The low frequencies are the ones that make the 
largest contribution to n. This is because high frequencies already have very short mean 
free paths, and so most of the heat in the alloy is carried by low frequency phonons. Thus, 
even though Si displays higher scattering rates at some intermediate and higher frequencies, 
it is the small scattering rates at low frequencies that determine the difference between the 
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FIG. 4: Thermal conductivity vs. nanoparticle diameter at 1% (top), 3.4% (middle), and 5% 
(bottom) nanoparticle concentrations, for temperature 800K. Triangles: Born+geometrical inter- 
polation. Squares: T matrix calculation for Si nanoparticle. Circles: T matrix calculation for Ge 
nanoparticle. 

lighter and heavier types of nanocomposite. 

A qualitative difference in the scattering rates of lighter and heavier impurities had been 
shown as early as 1963 for a model FCC scalar lattice with single and double impurities^i. 
Our results for the SiGe problem display the same kind of behavior, where the heavy impu- 
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FIG. 5: The ratio ngi/ nce, of the thermal conductivity of a SiGe matrix with embedded Si nanopar- 
ticles, Ksi, to the thermal conductivity of a SiGe matrix with embedded Ge nanoparticles, HQe, is 
shown as a function of nanoparticle diameter, at 300K, for three different nanoparticle concentra- 
tions, fp. 

rities scatter more strongly than the light ones at low frequency. This can be qualitatively 
understood by making an analogy with electron scattering by a local potential. Heavier 
impurities are analogous to a potential well, whereas lighter ones are analogous to a poten- 
tial barrier. From elementary scattering theory, the low frequency scattering cross section 
of a potential well is larger than that of a potential hump, in agreement with the trend 
observed. Nevertheless, to our knowledge, the effect of arbitrarily sized nanoparticles on a 
realistic 3 dimensional system has not previously been quantitatively investigated, and its 
consequences on thermal conductivity have not been addressed. 

Some further comments are in order. In principle, the techniques presented here would 
also allow us to study other nanoparticle shapes and compositions. Our choice of pure Si 
and Ge spherical nanoparticles has been motivated by simplicity. Experimentally, it may 
prove difficult to embed such nanoparticles into a SiGe matrix, because Si and Ge are fully 
miscible. Although high concentration Ge nanoparticles with fiat pyramidal or hemispherical 
shapes have been grown inside a Si^ and SiGe matrix— in the past, for the sake of clarity 
we have avoided introducing any experimentally determined morphological characteristics 
in our calculation. 
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IV. CONCLUSIONS 



We have preformed a parameter free first principles calculation of the thermal conductiv- 
ity of SiGe alloys with embedded Si or Ge nanoparticles. In contrast with the commonly used 
Born approximation, it is found that embedding nanoparticles in the material affects its ther- 
mal conductivity differently depending on whether the nanoparticles are relatively heavier or 
lighter than the embedding matrix. The calculation predicts that heavier nanoparticles (Ge) 
should be more efficient than lighter ones (Si) in reducing the k of Sio.sGeo.s- This behavior 
is determined by the higher scattering rate for heavier nanoparticles at low frequency, which 
is not predicted by the standard Born approximation, but is captured by the full Green's 
function calculation. Nevertheless, the approximated Born + geometrical approximation is 
found to work remarkably well, being within 20% of the exact result. The ab initio calcu- 
lation also confirms the existence of an optimal nanoparticle size that minimizes thermal 
conductivity, which had been previously predicted using a simpler model^. 
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